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SUMMARY 




Within the Born approximation and employing parabolic coordinates 
an expression is found which gives the cross section, induced by elec- 
tron collision, for excitation of the hydrogen atom from any initial state 
to any final state. Using this expression the cross section in the energy 
range of interest in plasma calculations for the following transitions 
among the principal quantum numbers are tabulated: n = 1, n -2,3,4, 

5, 6, 7, 8, 9, 10; n = 2, n' = 3, 4, 5, 6, 7, 8; n = 3, n' = 4, 5, 6, 7, 8; n - 4, 
n' = 5 , 6; n = 5 , n' = 6 . In conclusion a curve for the total inelastic col- 
lision of electrons and the hydrogen atom in its first five energy levels 


is constructed. 


I. INTRODUCTION 


The excitation cross section in hydrogen induced by electron col- 
lision calculated in the Born approximation is proportional to the 
squared modulus of the atomic form factor given by 

V(i, f) = je* 1 «/>. (r)0* (r)d 3 r , 

where and \fi f are the initial and final eigenfunctions of the atomic 
electron and K is the magnitude of momentum transfer of the incident 
electron. In this paper a closed form is found for the above expression 
when \p i and \p { are hydrogenic functions expressed in parabolic coor- 
dinates. Elwertl (1955) has evaluated this expression with similar 
specifications, although his final result is in differential form. 

The main concern of this paper is the evaluation of the electron 
impact induced excitation cross section between two arbitrary levels 
of hydrogen, calculated in the Born approximation. Up to now many 
calculations in the Born approximation have been carried out in this 
respect, and tables of cross sections with initial states in the range 
1-5 principal quantum numbers and final states 2-10 principal quantum 
numbers are available, 2 ' 10 although for higher levels the calculations 
are only for certain substates. 

In this paper, after formulation of the problem, the results in 
parabolic coordinates are compared with those in spherical coordinates 
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and their consistencies are examined. The calculation is then extended 


to higher levels, for which results are not available. All cross sections 
are listed in different tables. It is hoped that these tables will be 
useful in plasma and astrophysical calculations. 
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H. FORMULATION 


Excitation Amplitude 


Let the propagation vector of the exciting electron before and after 
collision be designated by k Q and kj, and the states of the atom in 
parabolic coordinates before and after collision by n 1 n 2 m and n x ' n 2 ' m. 
The excitation cross section in atomic units for such a collision is then 
given by 11 


/ § i i\ _ 8 tt f ° 1 i / , # x | 2 dK 

Q( n i n 2 m, nj n 2 m ) - — 5 | V(n l n 2 m, n x n 2 m) | — 

0 Jk <f k i 


(1) 


V ( n l n 2 m - n l' n 2 ">') ~ J eUCZ 4 (£ + 7 7> d V d 4> 


= s 


n °0 

3 ex P[T (f-’i)-'l (« + «') (f + ’)>] 


* (ft>* L ;,- (“f > L „>. <"?) L „*■*. (<*' ») (f + ij) df *, , (2) 


N n n being the normalization factor of the 77 eigenfunctions given 
by Bethe and Salpeter ,! 2 and also found in I. Similarly, ^ ni # n 2 # the 
factor corresponding to rj\ With this equation and the generating 
function of the associated Laguerre polynomials (cf. Eq. (37), I) it 
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follows that 
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where we have introduced 


P ” ~2~ ( a + a * ) 



(4) 
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a * s '\ 
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Before we carry out the differentiation of U with respect to p, we 
expand U in powers of s, s', t, t'. .Consider the expansion 


y(s, s') - (p + q + -rh + f~r) 0)s^s'^, 


^,= 0 ‘l 1=0 

(6) 

ix: 

y (0, 0) representing the and the derivatives of y(s, t) with 

't Z ' 

respect to s and t, evaluated at s = t n 0. To evaluate y 1 1 (s, s') , 
it is necessary to introduce 


as 


u 


1 - s ’ 


a’ s' 


(7) 


then 


y(u, v) - (p + q + u+v) ( m+1 ) 


( 8 ) 


It is convenient to introduce also 


g - (1 - s) , 


(1-s T 1 • 


(9) 


Then, making note of the relations 


du 

dl - a e' 


dv 


7 - a' h 2 , 


ds 


di~g 


n = ng n+1 


ds' 


h n = nh 


n+1 
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it follows that 


y 10 ( s , s') = ag 2 y' (u, v) , 

y 20 (s, s' ) = a 2 g 4 y 2 (u, v) + 2ag 3 y' (u, v) , 

y 30 ( s , s') = a 3 g 6 y 3 (u, v)+6a 2 g 5 y 2 (u,v)+ 6ag 4 y ' (u, v) , 

y v (u, v) being the v th derivative of y(u, v) with respect to either u or 
v. Inspection of the above equation shows that in general we can write 

l 

y^° (s, s') = C(v, 1) a v g^ +v y v (u, v) , 

V=1 

with C(v, -t) some constants. This equation is identical to Eq. (59), I, 
provided we let t — s' and a 2 - a in the latter equation. The dependence 
of y^° (s, s') on y v (u, v) is then similar to the dependence of y^(s, t) 
on y v (u, v) in I, and the C(y, satisfies the following recursion 
formula (cf. Eq. (60), I) 

C(v,^+1) = ('t + v)C(y ) 't)+C(y-l,'t), (10) 

with the boundary conditions , 

C(v, -£) - 0 when v = 0 or v>'i 1 , C(0, 0) = C(l, 1) = 1 . (H) 

A table of values of C(y, -C) is given in I. 
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Finally, corresponding to Eq. (61), I, or by direct deduction, we 


get 

't # 

Ml,. = V T ri„ 


(s, s') = 2_, 4 c h.*i) c K, 


v x =0 v x = 0 


( 12 ) 


Evaluation of y Vl 1 (u, v) when s - s' - 0 and substitution of Eq. (12) 
in Eq. (6) gives 


. ,\ "( m+1 ) 
as as 


l P + q+ !-• + 1- s' 



(-t 


. v,«; (n + v, +»,*)! 

, !<.' •) (-) ' » ' ini 


ViVi' 


x C(„, -t.) Cf,,' -t,') a-' a'"' X x s' 1 ' s'*' , (13) 


where we have introduced 


a - p + q 


(14) 


Similarly, 


. , . ,\ _ ( m+1 ) 
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(V* 2 ' !) (-) 2 2 x 


i. . v,+v,- (m+‘V 2 + v 2 ;- 


mT 


•V 


^2^ 2 V 2 V 2 
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With substitution of Eqs. (13, 15) in Eq. (5) we find 


au _ d_ d 

dp ~ \da + da*) ^ 


I 


‘1 . v 1 +v/+v 2 +v 2 # 


A '■ a ! t , ! a !) (-) 1 1 ! 2 x (» + ♦ A) ! (» + + a)’ 




V,+ V„ ,v[ +V 2 


x c(zW) C (^2 ^ 2 ) C (^ 2 ' ^2') a 1 2a ' 1 


x a 


-( m +2*' 1 *l'j) -( m +2 + v 2 +v 2 '> 


A A A „A' 


^|m + 1 + Vj + VjJ a*+|m+l + i' 2 + v 2 )ajs 1 s* 1 t 2 t' 

(16) 

The right-hand side of Eq. (3) becomes now, after making a bi- 
nomial expansion of the denominator, 


IN N , , 

1 n j n 2 n 1 112 
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X s 


H + ^l s ' j l + ^l t j 2 + ^2 t - h + ^2 
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Equating the coefficients of equal powers of s, s', t, t' of this equation 
and the left-hand side of Eq. (3), and substituting the value (cf. Eq. (20), I) 


N 


* 1“2 


1/2 


a 1 


ra+3/2 


n i ! n 2 ! 


+m) ! 3 (n 2 + m) 


i 3 


1/2 


a - 


n 


(17) 


we obtain 


V n 2 m, n 2 * n 2 ' mj 


72m+3 


n 1 ! n 2 ! ! n 2 ' ! 


1/2 


! (n 2 + m j ! ^rij' + m) ! (n 2 ' + m) ! 


L(T)(T) (^‘K^Kv vv.v.) <->— 




X .(m + Vj+v/)! jm + Vj+Vj'jlC^ych-tJC^'y C(v 2 ' 


a V '* V * a' V '* V * x a _ ( m+ 2 +v i +I/ i' ) x a *-(- + ^ 2 W). x ^ + l + Vj + a * 


+ + 1 + v 2 + a . (18) 


Here 7 stands for the set of 12 variable integers, 
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(j 1 j 2 i 1 i 2 ^l V 1 ^ 2 V 2 ^l V 1 ^2 V 2 ) 


(19) 
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subject to the restrictions 


o 

to 

i » 

ii 

- n j-'tj ; 

V 1 

= o, 1 , 2 , ••• , 


0, 1, 2, ••• 

’ * 

h 

= n 2 -^ 2 ; 

V 2 

= 0, 1, 2, • • • , 

*2 

0, 1, 2, ••• 

, ; 

h 



= 0, 1, 2, ••• , 


0, 1, 2, 


ia 

= n 2* ~^2 ; 

V 2 

= 0, 1, 2, , 





>( 20 ) 


in the transition n x n 2 ' m. Eq. (18) gives the desired excitation 

amplitude. 


Integration With Respect to K 


When Eq. (18) is substituted in Eq. (1), and the integration is car- 
ried out numerically with respect to K, the excitation cross section is 
obtained. In some cases it is advantageous to carry this integration 
analitically. To do this, we introduce 


/? - a + a ' ; 


( 21 ) 


then by Eqs. (4, 14), a = (/8 - iK) . We introduce further: 


r 


A = 2 (2 z ) 2m+3 (nn') 


-(m+2) 


n x ! n 2 ! ! n 2 ! 


1/2 


+ m) ! (n 2 + ! fn^ + ! (n 2 ' + m 


( 22 ) 
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G(y) = (-2) 


Vj+V 2 +V / +v 2 ' 


+ j A /m + j A /m + j A /m + j 2 ' ' 


x ! 4 2 ! V ! * a ' !) - 1 (» + v, + v/) ! (m + v 2 + v 2 ') ! C(v, 4J C(* 2 l 2 ) C(v/ V) C(* 2 ‘ * 2 ') 


x a*' 1 *' 2 a ,Vl+ " 2 ' (23) 


-(m+2+v +v /) 

H(y) = (/3-ik) 1 x ' 


„ -(m+2+v,+v ') 

x (/3 + iK) V 2 2 


x £(2m + 2 + Vj + v 2 + v i + v 2 'j /? + i - v 2 ~ v 2 'J kJ . (24) 

Then by Eqs. (1, 18) the cross section will be given by 


Q(' 
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8ttA 2 


G W G W IpWa). . 
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where N is the number of combinations in the set given by Eq. (19). 
By writing 


H (7) = (/32+K 2 ) 


■(m+2+Vj+i l 2 + i'| l +Vj ) 


x (/? + iK) 


v i* v i 


(/^ - iK) 


V 2 +V 2 
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1 

L 
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a(t)(iK)‘ 


= +K 2 ) -("' + 2 +v 1 + v 2 + Vi' + v 2 ') 
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where 


a(0) = (2m + 2 + v t +v 2 +i/ x ' + v 2 ') /3, a(l) - v i “ v 2 + ~ > ( 27 ) 

we get 


I (n < ? 2 ) = 2ZZ] ( _ > CT2(i)Vcr2 L ( y i’ "1) L ( 7 2 ’ "2) 


a. + cr 

K 1 2 dK 


K^+K*) 2 ”' 4 *^ ’ 
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where 
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(-) q 
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A.-p- 


p-q 
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A. - v i +v 2 +v i' +v 2 ’ a ~ p + q + 1 ; & - (p.q.t) , (29) 
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and each summation extends over all possible values of p, q, t . 
Noticing that ,1(72’ ^l) = I* (ti* ^ 2 ) * we can wr i te 


4ttA 2 


N N 


Q(n, n 2 m, n,' n,' m) - C(y,) G(y 2 ) [l(r,. > 2 ) * I*(),. 7 2 )] 

° (30) 


y i =1 r 2 = 1 


By Eq. (28), 


I (y 1 , y 2 ) + I *(y 1 , y 2 ) - [(~f 1 + C-)^] (i)^ 2 L( yi , L(y a , ^ 2 ) 


C7j +cr 2 


Jl 2 ~ 2, 2m + 4 + + \ 2 , / 3 2 , x x x 2 , (31) 


where we have defined 


i X 2 

x m (a + x)" n 


dx 


(32) 


Notice that expression (31) is always real. 


Symmetry Considerations 
It is evident from Eq. (2) that 

(33) 

(34) 


V ^n j n 2 ^ i n j n 2 ^ | -K^ - V*^n x n 2 m, n 2 ' m|K^ , 


V ^n 2 rij m, n 2 ' rij' m | ~K^ - n 2 m, r» 2 ' m |ltj 
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It follows that 


| V (rij n 2 m, n/ n 2 ' m|-K^ j 2 - | v ( n 2 n i m ’ n 2 # n i m | K ) | 2 

and, by Eq. (1), 


V( n i n 2 m, n x ' n 2 ' m 
(35) 


Q(n 2 m, n 2 ' m) - Q( n i n 2 n i # n 2 # m ) 


(36) 


Eqs. (35, 36) are used to test the accuracy of the numerical results. 


Multiplicity of States and the Total Cross Section 

Since the direction of the z-axis is taken along the momentum 
transfer vector K, the magnetic quantum number does not change in any 
transition. As n 1 +n 2 =n-m-l, n 1 can take the values 0,1,2,. . . , 
n-m-l;orn-m values. The same is true of n 2 . Then the total num- 
ber of combinations of and n 2 for a given n and m is n - m. Similarly, 
the total number of combinations of n^ and n' 2 for a given n' and m' is 


Designating the cross section for the transition nn^m-^n' n' 2 m 
by ^(nnj^m, n' n^n^m) , the cross section for the transition nn^m-^n'm 
is obtained by summing the former cross section over all the final states 
with a fixed m , 


n - m - 1 


(^(nnj^m, n' m) = ^ Q(nrij n 2 

n J = 0 


m. 


nn J n 2 


m), 


(37) 
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The cross section for the transition nm-*n'm is obtained by averaging 
(^(nnj^m, n'm) over all the initial states with a fixed m , 

n-m-l 

Q(nm, n'm) = (n - m) -1 ) Q(nn 1 r» 2 m, n'm). 

rijSsO 


(38) 


The cross section for the transition n -* n' is obtained by averaging 
Q(nm, n'm) with respect to the magnitude of the magnetic quantum number 


m. 


(39) 


n-1 

Q(n, n') = (2n - l)* 1 ^ ’ ^2 - 8(m,o)J Q(nm, n'm). 

m=0 

Since the total number of the initial states is 

n— 1 

j^2 - 8 (m, o)j (n - m) = n 2 , 

m=0 

Eq. (39) can be written alternatively as 

n~l n-m- l 

Q(n, n') = n -2 |2-8(m, o)J Q^n^m, n'm) # (41) 


(40) 


It is interesting to note that the number of independent transitions 
between the levels n and n 1 is given by 


n-l 


N = ) [ 2 - 8 (m, o) ] (n - m) (n' - m) . 

m=0 


When the right hand side is evaluated we obtain 


-tK- 


(42) 
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III. RESULTS AND DISCUSSION 


We have calculated within the Born approximation the excitation 
cross section of the hydrogen atom by electron collision for the tran- 
sitions n r 1 ton' = 2, 3, 4, 5, 6, 7, 8, 9, 10; n = 2 to n' = 3, 4, 5, 6, 7, 8; 

n - 3 to n ' “4,5,6,7,8;n = 4 ton' = 5, 6; and n = 5 to n' - 6 by 

employing parabolic coordinates. Previous similar calculations in the 
Born approximation using spherical coordinates have been made for 
the transitions n = 1 to n' - 2, 3, 4, 5, 6 by Me Carroll 4 ; = 2s to 
n' = 3, 4, 5, 6, 7 , 8, 9, 10 by Boyd 5 ; n4 = 2p , m = 0 , 1 to n 1 = 3, 4, 5, 6, 

7, 8, 9, 10 by McCrea and McKirgan 6 ; n - 3 to n' = 4 by McCoyd, 

Milford and Wahl 7 . There are few other calculations for certain op- 
tically allowed transitions between sublevels of higher levels but they 
do not give the total transition cross section between two levels. 


As a check on the consistency of the calculations, comparisons are 
made in this paper with all the values available in spherical coordinates 
Table I gives the excitation cross section from the ground state in both 
coordinates. The agreement is excellent. The cross section due to all 
higher states which are not listed explicitly can be calculated by a 


method given in Ref. 4. This is designated by 


CD 

L 

i=n+l 


Q(l, i ) > where n 


is the upper state of the highest transition whose cross section is 
listed in the table. Q(T) is the total excitation cross section. 


Table II compares the cross sections in the two coordinates for 
the transition n = 2 to n ' - 3 . Theoretically we must have 
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Q(2p ± 1) 


Q( 00, 01) + Q(00, 10) 


2Q( 00, 01) 


(43) 


Q(2s) + Q(2p 0) = Q(01, 02) + Q(01, 11) + Q(01, 20) 

+ Q(10, 02) + Q( 10, 11) + Q( 10, 20) 

= [2 Q(01, 02) + Q(01, 11) + Q( 01 , 20)] , (44) 

where Q(njn 2 , n t ' n 2 ' ) is the cross section for the transition between 
the sublevels n 1 n 2 and n 2 # n 2 ' , and where use has been made of the 
symmetry relation (36). The above equations are shown to be satisfied 
in Table II. Q(2, 3) is the excitation cross section for the transition 
n = 2 to n' = 3 , averaged over the initial states and summed over the 
final states. 

For brevity in the remaining of the tables only the averaged cross 
section for transition between the principal quantum numbers are 
listed. Table III gives the cross section for the transitions n ' 2 to 
n # - 4, 5, 6, 7, 8 , and the total excitation cross section from the n = 2 
level. Table IV is constructed to verify the excitation cross sections 
for n = 3 to n' = 4 as obtained in Ref. 7. Although the agreement is 
satisfactory, at low energy the two results differ to some extent. As 
our values are obtained both by the closed and the integral forms, and 
in the case of the integral form convergence to the numbers given is 
achieved by decreasing the mesh sizes of the numerical integration, 
it is believed that our results are more accurate. 

Table V gives the excitation cross section for the n - 3 to the 
n' = 4, 5, 6, 7, 8 levels, and the total excitation cross section from the 
n = 3 level. The contribution to the total cross section of the states 
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not listed is obtained by the method outlined in Ref. 4 and the known 
value of ionization of the n - 3 as given in I. 


To test the accuracy of the Born approximation, it is necessary to 
compare the result of the Born calculation with experiment. This is 
done in Fig. 1 along with the more elaborate theoretical calculation of 
close coupling 13,14 , and the classical theory of excitation given by 
Gryzinski 15 . According to this classical theory, if the energy of the 
incident electron is given as E 0 in rydberg and the atom is excited 
from the state n to n 1 , the excitation cross section is given by 


Q(n, n' ) - a n 4 a( 1 + a) 



g(n, n' + 1) 

[-fh )']' 


g(n, n') 


- ay 


+ (1 + a)y 


1 

3 ’ 


[- ay 2 + (1 + a)y] 3/2 


a 0 n 4.0307Tra 0 2 , a = (n 2 E 0 )- 1 , y 


ay £ 1 , 


ay > 1 , 

2 ' {$} 2 ■ <«> 


Compared to experiment, the Born approximation gives too high 
values, the classical theory gives too low, and the close coupling ap- 
proximation gives the best agreement. 

The disagreement between Born calculations and experiment may 
get worse for the excitation of higher states. This is due to the form 
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of the wave function of the bound electron. Since hydrogenic wave func- 
tion is used to evaluate the matrix elements of the Born calculations, it 
is implicitly assumed that the interaction potential between the two 
electrons is small compared with the interaction of the nucleus and 
the atomic electron. This, however, may not be the case for the excited 
states where the average distance of the electron from the nucleus 
is large. 

Fig. 2 compares different excitation cross sections and the ioniza- 
tion cross section of the n = 1 level. Figs. 3 and 4 make the same 
comparison for the n = 2 and the n n 3 levels. 

Table VI gives the excitation cross section of n = 4 ton' “ 5, 6; 
and the total excitation cross section for this level, while Table VII 
gives the excitation cross section of n = 5 to n' = 6 . Figs. 5 and 6 
show the results of these tables graphically. In Fig. 6 the classical 
curve is also drawn for comparison. The classical description of the 
excitation is open to question as transition to discrete levels cannot 
be described classically. 

Although the excitation of the ground state of hydrogen to the 2p state 
has an order of magnitude larger cross section than the excitation to 

I Q T 4 

the 2s state, especially at high incident electron energies it should 

be argued that this does not mean that for the excitation of higher 
states the cross section of non-optically allowed transitions can be 
neglected. This neglect is valid when Ka « 1* where K is the magni- 
tude of the momentum transfer of the incident electron and a is the 
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extent of the charge distribution of the bound electron before and after 
collision. For the excitation of high levels, a becomes more than an 
order of magnitude larger than the corresponding a for the ground state. 
This fact is evidenced by noting that the excitation cross section for 
optically allowed transitions in n = 5 to n' = 6 given in Ref. 9 is con- 
siderably lower than that given in Fig. 6. Similarly it can be argued 
that the dipole approximation within the Born approximation becomes 
less valid for excitation of the higher levels. 

Fig. 7 gives the total inelastic cross section by electron collision, 
including excitations to all levels and ionizations, for the first 5 levels 
of atomic hydrogen. 

In closing it should be mentioned that the cross section for de- 
excitation induced by electron collision is obtained through 

k i 

Q(f, i) = Q( i > f) . 

where Q(i, f) is the corresponding excitation, and k. and k f are the 
initial and final wave numbers of the incident electron in the excitation 
process. 

I wish to thank Robert Baxter for the efficient programming of this 
problem, Bernard Rugg for taking control over the production of the 
tables, and the Computing Branch of the Theoretical Division for pro- 
viding their generous services. 
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Table I. Excitation cross section of n = 1 level in units of rr . S is for spherical , and P 
is for parabolic, coordinates used in the calculation of the cross sections. 
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See Ref. 4. 

By interpolation. 



Table II. Excitation cross sections among the sublevels of n = 2 — 
n 1 = 3 transition in units of jt a^. m is the absolute value of the magnetic 
quantum number. Cross sections in parabolic coordinates are desig- 
nated by Q{n 1 n 2 , n' n^ ) with n x n 2 and n'n' belonging to the initial and 
the final states. 2^ is the sum of Q(2s) and Q(2p0), while 2^ is the sum 
of the three cross sections preceding^. Note has been taken of the 
relation Q(nj n 2 , n^ n' ) = Q(n 2 n x , n' nj ). The last two columns show 
the excitation cross section of n = 2 to n* = 3 in spherical and parabolic 
coordinates. 
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a. See Ref. 6. 

b. See Ref. 5. 




Table III. Excitation cross sections of n = 2 level to n' = 4, 5, 6, 7, 
8 levels in spherical and parabolic coordinates in units of tt a*. For 
spherical coordinates see Ref. 5,6. 





































































Table IV. Comparison of the excitation cross section for the tran- 
sition n = 3 to n 1 = 4 in spherical and parabolic coordinates in units of 
7 raj*. For spherical coordinates see Ref. 7. 
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Table V. Excitation cross section for the transition n = 3 to n 1 = 
4, 5, 6, 7, 8, in units of 77 a^. 




































































Table VI. Excitation cross section for the transition n = 4 to n' 


5, 6. 
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Table VII. Excitation cross section for the transition n = 5 to n* = 
6 in units of 77 a £ . 
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Fig. 2. Excitation of the ground state of the hydrogen to the n = 2, 3, 4, 
5, 6, 7 states. Q(l,i) is the ionization cross section of the ground state 
by electron collision. 
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